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ABSTRACT 

Recent observations of pulsar wind nebulae and radio polarization profiles 



O 



revealed a tendency of the alignment between the spin and velocity directions 
in neutron stars. We study the condition for spin-kick alignment using a toy 
g \ model, in which the kick consists of many off-centered, randomly-oriented thrusts. 

Both analytical considerations and numerical simulations indicate that spin-kick 
C\| \ alignment cannot be easily achieved if the proto-neutron star does not possess 

some initial angular momentum, contrary to some previous claims. To obtain 
the observed spin-kick misalignment angle distribution, the initial spin period of 
\ the neutron star must be smaller than the kick timescale. Typically, an initial 

^ \ period of a hundred milliseconds or less is required. 
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1. Introduction 

It is well known that pulsars have much larger space velocities than their progenitors, 
implying a kick at neutron star (NS) birth (e.g., Lorimer et al. 1997; Arzoumanian et al. 2002; 
Chatterjee et al. 2005; Hobbs et al. 2005; Winkler & Petre 2006). The physical mechanism 
for the kick, however, remains unclear (e.g., Lai 2004; Janka et al. 2005). One of the reasons 
that it has been difficult to constrain various kick mechanisms is the lack of correlations 
between kick velocity and the other properties of NSs. This situation has been changed 
due to the recent high- resolution Chandra X-ray observations of pulsar wind nebulae (e.g., 
Pavlov et al. 2000; Helfand et al. 2001; Ng & Romani 2004), which provided the evidence 
for spin-kick alignment for several young pulsars (e.g, Lai et al. 2001; Romani 2004; Wang 
et al. 2006) 1 . 
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1 A recent re-analysis of the proper motion of the Crab pulsar (Ng & Romani 2006) indicates that the 
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Another well-known method to determine the spin axis of pulsars is by the linear po- 
larization profile of radio emission. If the polarization profile could be described by rotating 
vector model (RVM), one can constrain the projected spin axis by the polarization angle at 
the center of the pulse. Previous attempts using this method have yielded ambiguous results 
(e.g., Morris et al. 1979; Anderson & Lyne 1983; Deshpande et al. 1999) mainly because 
for many pulsars the polarization profiles are not well described by the RVM (Weisberg et 
al. 1999). With Parkes surveys (e.g., Manchester et al. 1996, 2001), many well-calibrated 
polarization profiles became available, and more pulsar rotation measures to eliminate the 
Faraday rotation effect (e.g. Han et al. 2006). Moreover, proper motions for more than 200 
pulsars have been determined (Hobbs et al. 2005). By selecting pulsars with well calibrated 
polarization and proper motion measurements, Wang et al. (2006) have obtained spin-kick 
misalignment angle for 24 pulsars, and the data revealed a strong tendency of spin-kick align- 
ment. Johnston et al. (2005) independently obtained similar results for 25 pulsars based on 
different sample. 

On the other hand, one can constrain NS kicks using the orbital properties of NS 
binary systems (e.g., Dewey & Cordes 1987; Fryer & Kalogera 1997; Willem et al. 2004; 
Thorsett et al. 2005). In Wang et al. (2006), we obtained constraints on the kick magnitudes 
and directions for various NS binaries, including double NS systems, binaries with massive 
main-sequence companions, and binaries with massive white-dwarf companions. We found 
that the kick velocity is misaligned with the NS spin axis in a number of systems, and the 
NS spin period (when available) in these systems is generally longer than several hundreds 
milliseconds. 

What is the implications of the apparent spin-kick alignment for many pulsars? One 
possibility which is widely discussed (e.g., Johnston et al. 2005) comes from Spruit & Phinney 
(1998). They suggested that the initial spin of NS may originate from off-centered kicks even 
when the proto-NS has no angular momentum. They further suggested that if one imagine 
that the kick is composed of many random thrusts, then with multiple thrusts, alignment 
may be easily achieved. 

In this paper, we systematically study the condition of spin-kick alignment using a toy 
model similar to that of Spruit & Phinney (1998). We consider both the cases of zero and 
finite initial proto-NS spin. In §2 we introduce our toy model for kicks and describe our 
simulation procedure. In §3 we derive approximate but analytic conditions for spin-kick 
alignment. In §4, we present our simulation results (especially the distribution of spin-kick 
misalignment angle) under different initial conditions. We find that, consistently with our 
analytical estimate, without initial spin of the proto-NS, it is difficult to achieve spin-kick 
alignment. This is contrary to some previous claims (e.g., Spruit & Phinney 1998; Johnston 
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et al. 2005). However, with sufficiently short initial spin period (less than the timescale for 
each kick thrust), spin- kick alignment can be achieved. We discuss the implications of our 
results in §5. 



2. A Toy Model for Kicks 

The basic equations governing the evolution of the center-of-mass velocity V and angular 
velocity f2 of a proto-NS (mass M, radius Rns) are 

dV 

M^r = F, (1) 



dt 
dJ 

dt 



r x F. (2) 



Here F is the kick force, and r specifies the location where the force is applied, J = ICt is 
the angular momentum, with / = kMR^ s the moment of inertia. We adopt M = 1.4M , 
i? NS = 10 km and k = 0.4. 

We model the kick force on the NS as consisting of n thrusts, Fj, % — 1, 2, n, each has 
a duration Tj. During each thrust, we construct a "temporary" body frame (xyz) corotating 
with the star so that Fj is constant in this frame (see Fig. 1). We specify the kick force 
Fj by the magnitude Fi (= F, the same for all thrusts) and two angles aj,/3j. We choose 
(ctj, f3i) randomly distributed in the range of of 0° < «j < 30° and 0° < $ < 360°. Each 
thrust acts at the position rj, which is specified by the spherical coordinates (rj,0j,0j) in 
the nonrotating frame (XYZ) at the beginning of the thrust. Note that choosing a different 
range of ctj is equivalent to choosing different rj, as long as rjsinctj remains the same. 
During the course of each thrust (duration Tj), the rj changes in the XYZ frame as the body 
rotates. In our simulation, we assume Tj = r is the same for all thrusts. Thus the total 
kick time is Tk; c k = ^ r « = nT - We set the total momentum V = FT^^ at the fixed range 
V = MV = M(500 - 2000) km s" 1 in all simulations. We fix n = 30 km for all thrusts. 
For (#j,0j), we consider two possibilities: (i) The kick position is randomly distributed on a 
sphere in the XYZ frame; (ii) The kick position is randomly distributed on a sphere in the 
body frame of the NS. Our simulation results reported in §4 refer to the first case. We have 
found that the results of the second case are similar. 

For each thrust, r^ and Fj are fixed in the the "temporary" body frame, i.e., rj = TjZ, 
and Fj = -Fj(sin «j cos /3jX + sin «j sin /3jy + cos ttjz). To find the components of the force Fj 
and torque rj x Fj in the inertial frame (XYZ), we need to solve for the time evolution of the 
body axes e (= x, y or z). The (XYZ) frame and (xyz) frame are related by the rotation 
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Fig. 1. — Geometric model used in our simulations. XYZ is a nonrotating "fixed" frame 
centered at the neutron star. A thrust force F; is applied at the position r; (specified by the 
spherical coordinates r iy 9i, (pi. The body frame xyz is constructed with the z-axis along r,, 
and the x-axis in the meridional direction. The direction of Fj is specified by the two polar 
angles ctj, 
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matrix T: 











(I) 






l-( 












(3) 



So we have the expressions of the body axes in the (XYZ) frame: 






(4) 



At the beginning of each thrust 



cos 9j cos cf>i — sin <pi sin 9j cos <_ 
T= cos^sin0j cos0j sin 6^ sin 0j . (5) 
— sin $i cos 6i 




The body axis e evolves according to 



de 

dt 



SI x e. (6) 



Substituting Eq. (4) in Eq. (6), we obtain the evolution of each component of T, e.g., 
dTn/dt = Sl 2 T 31 — ^3^21, et al., where ^1,2,3 are the three components of fl in the XYZ 
frame. 

Using Eqs. (1) - (6), we can directly simulate the movement and the rotation of the NS. 
Consider a star with initial velocity V = V init and angular velocity Ho = ^init- Suppose it 
receives a thrust Fi at a random position ri with duration T\. We use the 4th order Runge- 
Kutta method to integrate equations (1), (2) and (6) to obtain Vi and 1; the velocity and 
rotation rate after the first thrust. If the star receives a new thrust, we just select a new 
body frame according to the new position of the thrust, considering Vi and fl\ as the initial 
velocity and angular velocity, and repeat the calculation as in the first thrust and so on. 



3. Analytic Consideration 

Consider a proto-NS without initial spin. The kick consists of n thrusts, each with the 
same duration r. We assume n ^> 1 in this section. After the first thrust, the star receives 
an angular velocity 

. ^ ^ Frr sin a Frr sin a 

AQ = = = 5—, (7) 

/ kMR 2 NS ' V ; 
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and a velocity AV not more than Ft/M. Note that after the first thrust, the star's spin fii 
is always perpendicular to its velocity V x . If the duration of each thrust is larger than the 
spin period caused by the first thrust, i.e. 

then the second thrust will be rotationally averaged such that the net thrust will be along 
f2i (the unit vector along fli). similar argument applies to additional thrusts. The final 
characteristic velocity and angular velocity are given by 

<Vf> = ((£ AV *-ni) )~^ 2 > (9) 

(£l 2 f ) = / A«i -fill \ ~ ^nAQ 2 . (10) 
The typical final velocity and angular velocity are aligned, with 

Vf^J^Avn,, Qf-J^Ann^ (n) 



Equation (8) is a sufficient condition for spin-kick alignment but not a necessary one. 
It's convenient to define a critical ratio 

^{M- (i2) 

For r < 2n/AQ or n c > 1, spin-kick alignment may or may not be achieved. If n c > n or 
r < 2n/ (y/nAQ), different thrusts add up in a random walk fashion. The final spin and kick 
are of order 

V f ~ VnAV, n f ~ v^A^, (13) 
with random angle between the spin and kick. 

For 1 n c < n, or 2ir / (y/nAVl) < r ^ 27r/Af2, the situation is more complicated. For 
the first n c thrusts (i = 1, 2, n c ), the thrust duration r satisfies r < 2ir / '(^/iAil) , thus 
the characteristic velocity and rotation rate are 

Vi-ViAV, Qi~ViAn, i = l,2,---,n c (14) 

with random directions between Vj and f2j. For the remaining thrusts (i — n c + 1, n c + 2, 
n), r > 27r/f2j, so that rotationally averaging is effective. The final velocity and angular 
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velocity are 



V/ ~ V nc + J^—^AVn nc , (15) 



n — n c - / n — n c 



n f ~ n„ c + y— ^Afifi nc ~ _£j AQ^ nc . (16) 

Here V„ c ~ ^fn~ c lSV has random direction compared to £l nc . So spin and kick will be aligned 
when \/ (n — n c )/3 ^> which means 

n > 4n c . (17) 

Otherwise, spin and kick will be misaligned. 

To summarize, spin-kick alignment /misalignment depends on the critical ratio n c (see 
Eq. 12). Let F = V/T kick = MV/T kick , r = f n R NS , we find 

2irkR NS n \ 2 ( n 2 , ns . 



_ fnrV sin a J \ 40/aVi fe Ti sin a J 
where we have used T kick = nr, T-y = T kick /(1 s), and V\k = V/(10 3 km s -1 ). For 

n c <max(n/4, 1), (19) 

spin and kick will be aligned, while for n c > max(n/4, 1), spin and kick will be misaligned. 

If the NS has initial spin H init , the sufficient alignment condition, equation (8), should 
be modified to 

2tt 

T ~ max(Afi, fi init )' ^ 



4. Simulation Results 

In our model, the key parameters are n and r or T kick , as well as the initial spin period 
Pinit- Depending on the kick mechanisms, the total kick duration T kick ranges from 0.1 s to 
a few seconds (e.g., Lai et al. 2001; Socrates et al. 2005; Scheck et al. 2006; Burrows et 
al. 2006a, b). Note that we choose V = FT kick in the range of M(500 - 2000)kms" 1 and 
other parameters such that the final distributions of kick velocity and spin period of NSs 
qualitatively agree with observations (see Hobbs et al. 2005). 



In Fig. 2 - 4, we present simulations of 20000 pulsars without initial spin (ilmit = 0). 
The final spins of the NSs are all due to the off-centered thrusts. In Fig. 2, we fix the total 
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kick duration to T kick = 1 s, while changing the number of thrusts: n — 5, 10, 20. For these 
cases, n c ~ [n/ (7y/Vik)\ (see Eq. 19). We find that an aligned distribution is produced for 
n — 5, but not for n = 10 or 20. 

In Fig 3, we consider different values of total kick duration T kick = 0.1 s, 0.5 s, 1 s, while 
fixing the number of thrusts to n — 5. For these cases, n c ~ .2 / (V^Ti) 2 . So we find that 
for large T kick , an aligned 7 distribution is produced. 

Figure 4 shows the cases with the same thrust duration r = 0.2 s, while the number 
of thrusts are n — 5, 10, 20. Here n c ~ [n/(60V lk Ti)] 2 (where T\ = r/1 s). Since either 
n c < 1 or n c < n/4 is satisfied for these cases, the 7 distributions all show an tendency of 
alignment. 

Note that in the above three figures, the kick velocities are all distributed at a few 
hundred kilometers per second and the final spin periods are distributed from 10 to hundreds 
of milliseconds, in agreement with observations. Although an aligned 7 distribution can be 
produced under certain conditions (see Eq. 19) without initial spin, the distributions are 
significantly broader than what was observed (Johnston et al. 2005; Wang et al. 2006; see 
Fig. 5). 

Figure 6 shows the simulations with different initial spin period Pi n i t = 500 ms, 100 ms 
and 50ms, all with T kick = Is and n = 10. Clearly, for P init < r, rotational averaging 
is effective, and spin-kick alignment is easily achieved. Figure 5 compares our simulation 
results with the observed spin-kick misalignment angles based on pulsar polarization profiles 
(see Wang et al. 2006), taking into account of the sky projection effect. We see that for the 
-Pinit = 50 ms simulation depicted in Fig. 6, the simulated spin-kick distribution agrees with 
observational data. With P init = 100ms (other parameters being the same), the simulated 
distribution is broader than the data. The key condition for producing alignment is Eq. (20). 

Figures 7 and 8 show the distribution of the misalignment angle 7 as a function of V) and 
Pf. We see that when the 7 distribution is broad (Fig. 7), pulsars with different V/'s have 
similar range of 7's. On the other hand, based on our toy model simulations, for an aligned 
7 distribution (Fig. 8), high-V/ pulsars have a strong tendency for spin-kick alignment. We 
have attempted to test such Vf — 7 correlation in the existing sample for 24 pulsars. The 
current data does not show such correlation, probably because of the small sample or large 
error in various measurements. 
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5. Conclusion 

In this paper, we have developed a toy model to study the conditions for pulsar spin- 
velocity alignment in supernova kicks. We have focused on the idea (Spruit & Phinney 
1998) that multiple off-centered thrusts to the proto-neutron star may result in spin-kick 
alignment. We found that without initial angular momentum, the alignment cannot be 
easily produced. To obtain the observed spin-kick alignment distribution based on radio 
pulsar polarization data (Johnston et al. 2005; Wang et al. 2006), the proton-neutron stars 
should have appreciable rotation rate, with period less than the timescale of each kick thrust. 
The typical initial period required is < 100 ms. 

Currently, the most conservative (and promising) kick mechanisms are "hydrodynami- 
cally driven kicks" . In particular, large-scale convections, instabilities or wave modes devel- 
oped in the neutrino-heated mantle behind the shock and in the proto-neutron star may natu- 
rally lead to asymmetric explosion (e.g., Thompson 2000; Scheck et al. 2004,2006; Blondin & 
Mezzacappa 2006; Foglizzo et al. 2005; Burrows et al. 2006a, b; Yamasaki & Yamada 2006). 
The supernova simulations cited above do not include initial angular momentum, and the 
resulting kicks are randomly distributed. It is possible that with even a small rotation, the 
hydrodynamical instabilities may preferentially develop along the rotation axis. If the kick 
timescale is long (as indicated by recent simulations), rotational averaging may be effective 
and a preferentially aligned spin-kick distribution can be produced. 

Other kick mechanisms (such as those based on asymmetric neutrino emissions in the 
proto-neutron star; e.g. Duncan & Thompson 1992; Lai & Qian 1998; Arras & Lai 1999a,b; 
Socrates et al. 2005) and the "electromagnetic rocket" effect (Harrison & Tademaru 1975; Lai 
et al. 2001) can easily result in spin-kick alignment, but they require more extreme conditions 
(such as superstrong magnetic field or very rapid spin) for the proto-neutron stars (see Lai 
2004; Wang et al. 2006 and references therein). 

This work is supported by National Natural Science Foundation of China (10328305, 
10473015 and 10521001). D.L. has also been supported in part by NSF grant AST 0307252 
and NASA grant NAG 5-12034. D.L. thanks NAOC (Beijing) for hospitality during the 
course of the work. 
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Fig. 2. — The distribution of the spin-kick misaligned angle 7 (upper panel), final velocity 
Vf (lower-left panel) and final spin period Pf (lower-right panel) of NSs in simulations with 
TLick = 1 s and different number of thrusts n =5, 10, 20 (from thick lines to thin lines). The 
initial spin of the NS is fl^t = 0. 
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Fig. 3. — Same as Fig. 2, except for different values of the total kick duration Tki C k =0.1s, 
0.5s, Is (from thick lines to thin lines). The number of thrusts is fixed at n — 5, and the 
initial spin is f2 init = 0. 
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Fig. 4. — Same as Fig. 2, except that the duration of each thrust is fixed at r = 0.2 s, and 
the number of thursts are n — 5, 10, 20 (from thick lines to thin lines). The other parameters 
keep the same as Fig. 2. 
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Projected Spin— Kick Angle 

Fig. 5. — The spin-kick alignment distributions from observations and simulaions. The 
points (with error bars) are based on radio pulsar polarization and proper motion data (see 
Fig. 1 in Wang et al. 2006, note that the spin-kick angle for the Crab pulsar has been changed 
from 8° ± 20° to 26° ± 3°, based on a recent analysis, see Ng & Romani 2006): each point 
represents a pulsar with measured angle between the projected spin axis and the proper 
motion direction. Note that because of the orthorgonal mode phenomena, there is a 90° 
degeneracy for the spin axis inferred from the pulsar polarization profile. Thus we have 
folded Fig. 1 of Wang et al. around 45° (e.g., a 60° data point is identified with 30°). The 
histograms are based on our simulations with P init = 50 ms, T kick = 1 s and n — 10 (see 
Fig. 6): The dashed line gives the actual (unprojected) spin-kick angle distribution, the solid 
and dot-dashed lines give the projected spin-kick angles in the plane of the sky when the 
line of sight is inclined with respect to the initial spin axis by 60° and 30°, respectively. To 
compare with observational data, the simuation results are also folded around 45°. 
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Fig. 6. — Same as Fig. 2, except for different values of the initial spin period Pi n it =50 ms, 
100 ms, 500 ms (from thick lines to thin lines). The other parameters are = 1 s, n = 10. 
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Fig. 7. — The distribution of the spin-kick misalignment angle 7 as a function of the final 
velocity Vf (left panel) and the period Pf (right panel) for a simulation which produces a 
broad 7 distribution. The simulation parameters are = 1 s, N = 20, fi; n i t = 0. 
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Fig. 8. — The distribution of the spin- kick misalignment angle 7 as a function of final 
velocity Vf (left panel) and period Pf (right panel) for a simulation which produces an 
aligned 7 distribution. The simulation parameters are — Is, N — 10, Pi n ; t = 50ms. 



